Synchronization from Disordered Driving Forces in Arrays of Coupled Oscillators 
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The effects of disorder in external forces on the dynamical behavior of coupled nonlinear oscillator 
networks are studied. When driven synchronously, i.e., all driving forces have the same phase, the 
networks display chaotic dynamics. We show that random phases in the driving forces result in 
regular, periodic network behavior. Intermediate phase disorder can produce network synchrony. 
Specifically, there is an optimal amount of phase disorder, which can induce the highest level of 
synchrony. These results demonstrate that the spatiotemporal structure of external influences can 
control chaos and lead to synchronization in nonlinear systems. 
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Networks of coupled nonlinear oscillators provide use- 
ful model systems for the study of a variety of phenomena 
in physics and biology |l|. Among many others, exam- 
ples from physics include solid state lasers Q and coupled 
Josephson junctions In biology, the central nervous 

system can be described as a complex network of oscil- 
lators [5j , and cultured networks of heart cells are exam- 
ples of biological structures with strong nearest-neighbor 
coupling |fg. In particular, the emergence of synchrony in 
such networks |7y8j and the control of chaos in nonlinear 
systems 0, TToL have received increased attention in 
recent years. 

Disorder and noise in physical systems usually tend 
to destroy spatial and temporal regularity. However, in 
nonlinear systems, often the opposite effect is found and 
intrinsically disordered processes, such as thermal fluc- 
tuations or mechanically randomized scattering, lead to 
surprisingly ordered patterns [lj. For instance, in the 
phenomenon of stochastic resonance the presence of noise 
can improve the ability of a system to transfer informa- 
tion reliably |l3j . Some time ago, Braiman et al. stud- 
ied one- (ID) and two-dimensional (2D) coupled arrays 
of forced, damped, nonlinear pendula |l4| . They found 
that when a certain amount of disorder was introduced 
by randomizing the lengths of the pendula the dynamics 
of the array ceased to be chaotic. Instead, they observed 
complex, yet regular, spatiotemporal patterns. Further 
studies of the same system showed that chaos in the ar- 
ray of oscillators can also be tamed by impurities El 
and that random shortcuts between the pendula lead to 
synchronization of the array |16| . 

Here, we introduce disorder by modifying the driving 
forces of the oscillators through phase differences. We 
observe the emergence of regular, phase-locked dynamics. 
Moreover, for intermediate spreads of the phase angles in 
the driving forces, we find that the oscillations become 
largely synchronous. 

We focus our numerical analysis on arrays of forced, 
damped, nonlinear pendula. The ID array (chain) is de- 
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scribed by the equation of motion 

ml 2 8 n + 7#„ = — mgl sin n + r' + r sin (ujt + ip n ) 

+K{6 n+1 + e n - 1 -2e n ), n = i,2,... n. (i) 

In order to consider a 2D lattice, we introduce an addi- 
tional index, n — > 9 n>m , n, m = 1, 2, . . . N and modify 
the coupling term accordingly: + 6>„_i — 20 n ) — > 

n(0 n +l,m + 8n-l,m + #n,m+l + #n,ro-l — ^n,m)- F° r both 

the ID and 2D case, we choose free boundary condi- 
tions, i.e., 6q — 6*i, 8n — On+i and o ,m = #i,m, On,™, = 
0N+i,m, Onfi = n ,i, ®n,N = n ,N+i, respectively. The 
parameter values used are the same as in previous stud- 
ies |3 EI El : The mass of the pendulum bob is m = 1 , 
the length 1 = 1, the acceleration due to gravity g = 1, 
the damping 7 = 0.75, the d.c. torque r' = 0.7155, the 
a.c. torque r = 0.4, the angular frequency to = 0.25, and 
the coupling strength k = 0.5. For this choice of pa- 
rameter values, each isolated pendulum displays chaotic 
behavior characterized by a positive Lyapunov exponent 
0. 

A particularly easy and intuitive way to visualize the 
global spatiotemporal behavior of a chain (or lattice) of 
oscillators is to consider the average velocity 

1 N 

^^E^w ( 2 ) 

at times that are integer multiples of the forcing period 
T = 2n/u> [15j. Considering this measure for an isolated 
pendulum, Gavrielides et al. perfomed a bifurcation anal- 
ysis with respect to the pendulum length I and found that 
an uncoupled pendulum is chaotic for values I — 1 ±0.002 
If the length of an isolated pendulum is increased to 
I > 1.002, it performs a 'libration,' in which the combined 
d.c. and a.c. torque are insufficient to overcome the pen- 
dulum's increased rotational inertia. On the other hand, 
if the pendulum's length is decreased to / < 0.998, the 
pendulum performs a 'rotation,' an overturning motion 
where the torques combine to rotate the pendulum over 
the top. 

In our study, we do not alter any parameters that 
would affect the dynamics of an isolated pendulum and 
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FIG. 2: Chaotic and regular dynamics as a function of the 
degree of disorder. The average angular velocity at t — 
60T, 6 IT, . . . , 80T is shown for each value of the disorder pa- 
rameter k. a) ID array of N — 50 oscillators, b) 2D lattice 
of 16 x 16 oscillators. 



FIG. 1: Spatiotemporal angular velocity plots for chaotic and 
regular dynamics in an array of N = 50 coupled oscillators. 
The chain of pendula is shown from left to right. Time in- 
creases continuously from bottom to top. Grayscales indicate 
the angular velocities of the oscillators. Light gray shades 
represent negative, dark tones positive velocities. 



keep the coupling strength at its default value. Instead, 
we introduce disorder by randomly varying the phase an- 
gles ip n of the driving forces in Eq. Q . In the case where 
(f n = for all driving forces, we observe chaotic dynamics 
in the array (Fig. ^) in agreement with previous studies 
|14|. However, when we disorder the driving forces by 
randomly choosing the phase angles ip n uniformly from 
the interval [— kir, -\-kn], we observe that for sufficiently 
large k the oscillations become regular. 

Figure [3] shows the average angular velocity cr(t) at 
t = 60T, 61T, . . . , 80T for a ID array of N = 50 and a 
2D lattice of 16 x 16 oscillators. The presence of chaos 
for small disorder in both the ID and 2D array becomes 
manifest in a dispersed distribution of the average veloc- 
ities cr(60T), cr(61T), . . . , cr(80T). For larger disorder, 
however, we observe periodic patterns in the form of 1T-, 
2T-, 3T-, . . . 'attractors,' where the average velocity of 
the oscillator array repeats its value after 1,2,3,... forc- 
ing periods. Ultimately, as k is increased further, a IT 
periodic pattern is reached. 

In general, the value of k for which a transition from 
chaotic to regular dynamics first occurs depends on the 



particular distribution of the random phases. We thus 
consider the average over several different samplings of 
uniform distributions in order to analyze the occurrence 
of different forms of periodic behavior. Figure [3] shows 
the probability for a ID array to have reached a IT-, 
2T-, 3T-, or 4T-attractor after t = 60T as a function of 
the disorder parameter k. For very small disorder, i.e., 
k < 0.02, we observe only chaotic dynamics, but as k 
passes this threshold, the first periodic patterns start to 
appear. For k > 0.1, we observe that IT-, 2T-, 3T-, 
AT-, . . . attractors coexist with chaotic behavior. For 
0.02 < k < 0.13 the 2T-attractor is the dominant form of 
dynamics if an attractor has been reached. For k > 0.28, 
the array undergoes regular oscillations with period IT 
in the vast majority of cases. 

Furthermore, in addition to the transition from chaotic 
to regular behavior, we observe that the oscillations be- 
come largely synchronous, i.e., the phases of the oscilla- 
tions not only lock but tend to assume equal values, for 
intermediate values of k. In order to quantify the pres- 
ence of synchrony in the array, we consider the averaged 
cross correlation 
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where Cij denotes the correlation between the ith and jth 
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FIG. 3: Probability P of chaotic dynamics (solid line) and dif- 
ferent forms of regular behavior (dashed lines) vs. the disorder 
parameter k in an array of N = 50 coupled oscillators. The 
probabilities were determined by averaging over 100 different 
samplings of the phases <p n . 
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FIG. 4: Synchronization in ID (dots) and 2D (squares) arrays 
of oscillators vs. the disorder parameter k. Error bars show 
one standard error of the mean. Filled dots correspond to 
N — 16, open dots to N = 50. Averaging was performed over 
200 (ID) and 10 (2D) different samplings of the phases ifi n - 



oscillator: 



S +T dtef(t)j^dt9](t) 



T Q +T 
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FigureQ]shows C as a function of k for two ID and one 2D 
arrays. Disordering the driving forces results in less syn- 
chronized oscillations of the array if the disorder param- 
eter is very small. The minimum of synchrony is reached 
for k w 0.03. Note that the location of this minimum cor- 
responds approximately to the first appearance of regular 
dynamics in Fig. [3J When the external forces are disor- 
dered further, synchronization in the array increases and 
reaches a peak value for intermediate disorder. In the ID 
case, the maximum is reached for k w 0.3 and its value is 
C* max ~ 0.72 for N = 50 and C max w 0.78 for N = 16 os- 
cillators. In the case of the 2D array, the synchronization 
is even stronger. Here, the peak value of C max ~ 0.95 
is reached for k » 0.2. We attribute the stronger syn- 
chronization in the 2D array to the fact that the num- 
ber of couplings per oscillator is higher than in the ID 
case. Furthermore, smaller arrays show a higher degree 
of averaged cross correlation than larger arrays. This is 
because oscillators that are nearest neighbors show the 
highest degree of synchronization, and the ratio of cross- 
correlation coefficients obtained from direct neighbors to 
all cross-correlation coefficients contributing to the aver- 
aged cross correlation C decreases with increasing size of 
the array like 0(l/N). 

To summarize, we have shown that disorder leads to 
transitions from chaotic to regular behavior in arrays of 
coupled oscillators when disorder is introduced in the 
phases of the driving forces ^^|. In this investigation, 
each pendulum was in a regime where it behaves chaoti- 



cally when uncoupled, in contrast to previous studies in 
which parameters were altered that affect the dynamics 
of an isolated oscillator |!4lll5l |. In particular, Braiman et 
al. introduced disorder by randomly varying the lengths 
of the pendula Since an isolated pendulum only 

behaves chaotically when its length lies within a narrow 
range, only 2% of the oscillators remained in their chaotic 
regime in this approach, and the transition from chaotic 
to regular spatiotemporal patterns reported in Ref. fl4| 
can be attributed to the dominance of the majority of 
regular pendula over the few remaining chaotic ones |15| . 
Our results show that disorder in the model system de- 
scribed by Eq. JQ) results in regular dynamics of the array 
even if all individual elements are chaotic. Moreover, we 
find that for intermediate disorder, the oscillations show 
a high degree of synchronization. 

Stimulus-induced synchronization of neural activity in 
central nervous systems has intrigued neuroscientists for 
decades ^| |2(j. Furthermore, in many applications, 
such as in coupled Josephson junctions, or in the case 
of atrial or ventricular fibrillation, one seeks to restore 
periodic or steady-state behavior from chaos. It is in re- 
gard to these day-to-day circumstances that control and 
synchronization of chaotic dynamics have become one of 
the central topics of nonlinear science [2ll |22^ | . In most 
situations the components of a system themselves can- 
not be altered, so it is desirable to establish methods by 
which chaos can be tamed without changing parameters 
intrinsic to the system. We thus believe that our pro- 
posed mechanism of controlling chaos via external forces 
has potential applications in these fields. 
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